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Operation  of  the  Facility  Model 

By:  C.  Haldeman  and  M.  Dunn 
Date:  5-7-93 

Introduction: 

The  facility  model  is  a  major  operating  tool  for  the  ATARR  facility.  It  can  provide 
information  about  setting  the  eddy-brake,  main  valve  activation  pressures,  initial  fill  pressures,  test 
gas  properties  as  well  as  predicting  performance.  It  is  important  to  remember  that  the  facility 
model  is  just  a  model,  the  numbers  it  produces  are  approximate.  Some  time  has  been  spent  with 
the  model  trying  to  make  the  predicted  main  valve  performance  reflect  its  actual  motion  more 
accurately.  In  the  future,  as  systems  are  added,  it  should  be  possible  to  compare  the  model 
predictions  of  the  torque,  eddy  brake  temperature,  etc.  with  data.  As  a  result,  the  facility  model 
will  generally  be  one  which  changes  with  time.  This  memorandum  is  written  to  document  some  of 
the  history  of  the  model,  discuss  how  it  is  run  presently,  and  point  towards  future  work. 

History: 

The  facility  model  (originally  called  the  math  model)  was  built  at  MIT  by  Dr.  Gerry 
Guenette.  The  model  arrived  at  Calspan  and  was  documented  initially  by  Mike  Stanton  (see 
included  "Original  Math  Model  Description"  by  Mike  Stanton),  and  not  much  was  done  with  it 
until  the  main  valve  model  verification  began.  Several  changes  were  made  in  the  default  values  for 
the  plumbing  sizes,  and  the  friction  coefficients  (to  reflect  reality).  These  are  all  documented  in 
three  memorandums'  attached  at  the  end  of  this  document.  All  of  these  changes  dealt  with  the  main 
valve  performance  aspects  of  the  model. 

Some  relatively  strange  behavior  has  been  observed  in  the  valve  performance.  This  has 
required  the  examination  of  the  original  main  valve  model  (which  was  a  four  chamber  model)  with 
a  view  towards  incorporating  two  more  choke  points.  Presently  WPAFB  is  performing  this  task. 
In  addition,  at  high  temperatures,  the  main  valve  did  not  operate  smoothly,  requiring  more  actuator 
pressure  to  open.  Upon  inspection,  it  was  determined  that  the  thermal  expansion  of  the  original 
bearings  was  incompatible  with  the  surrounding  material  and  that  some  binding  could  occur  at 
higher  temperatures.  These  bearings  have  been  redesigned  and  rebuilt,  and  testing  of  the  valve 
should  begin  shortly. 


'  C.  Haldeman  and  M.  Dunn,  "Main  Valve  Model  Update",  Memo  7-27-92;  C.  Haldeman  and  M.  Dunn,  "Main 
Valve  Model  Update  H",  Memo  10-6-92;  C.  Haldeman,  "Math  Model  Revisions",  Memo  10-27-92 
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In  all  other  respects,  the  facility  model  is  the  same  as  the  one  delivered  from  MU.  What 
has  been  changed  is  the  manner  in  which  the  model  is  run,  and  the  interpretation  of  the  results. 

For  detailed  discussion  about  the  integration  scheme,  and  the  governing  equations  of  the  model  one 
should  see  either  Dr.  Guenette's  notes  or  Mike  Stanton's  documentation.  The  rest  of  this 
document  will  discuss  how  the  model  is  operated  presently. 

Operation  of  Facility  Model 

The  entire  facility  model  suite  of  software  consists  of  a  main  FORTRAN  code  named 
"facility .F"  which  resides  in  /usr/atarr/src/fortran ,  collections  of  subroutines  called  by  "facility" 
(also  contained  in  the  same  directory),  and  an  X-1 1  interface  called  mm  (which  resides  in 
/usr/atarr/src/mathmodel).  Listed  in  "facility"  are  all  of  the  subroutines  which  the  main  code  calls. 

It  should  be  noted  that  all  executable  code  placed  in  /usr/atarr/bin  can  be  enacted  in  any  directory 
and  this  is  where  copies  of  all  DAS  system  executable  code  exist.  The  program  can  be  run  either 
as  a  stand-alone  program  by  typing  "facility"  or  using  an  X-1 1  interface  invoked  by  typing  "mm". 
In  both  cases  the  program  requires  at  least  two  files:  a  parameters  file  and  a  configuration  file.  The 
FORTRAN  code  presupposes  that  those  already  exist,  the  X-1 1  interface  allows  you  to  make  these 
files.  The  variables  which  are  listed  in  these  files  are  included  in  a  note  called  "Facility  Model 
Parameters",  located  at  the  back  of  this  report. 

Output  of  the  program  consists  of  several  files.  The  original  files  were  BD.OUT  and 
BD.PLT.  These  files  regurgitate  the  parameter  and  configuration  files  as  well  as  the  time  histories 
of  various  properties  (in  ASCII  format).  The  only  difference  in  these  files  is  that  the  BD.PLT  file 
was  designed  to  be  used  with  a  plotting  routine.  New  ouq)ut  files  have  been  added.  One  is 
defined  by  the  user  as  a  comma  separated  file  containing  main  valve  data  only  (default  name  is 
valve.dat).  This  is  a  consolidation  of  information  contained  in  BD.PLT  which  has  been  used  to 
analysis  the  main  valve  performance.  This  file  is  easily  portable  to  other  operating  systems  and 
other  programs  for  more  detailed  analysis.  In  addition,  there  are  several  files  of  the  form  P(5- 
8)_FM.3,  POS_FM.3,  and  VEL_FM.3.  These  are  time  history  files  of  the  main  valve  pressures 
and  motion  (position  and  velocity)  which  have  been  placed  into  a  DAS  context  (i.e.  they  conform 
to  the  DAS  architecture).  These  files  can  be  directly  read  by  DRP  (the  data  reduction  package)  and 
plotted  against  data  making  model  vs.  data  comparison  quite  easy. 

The  program  has  two  modes  of  operation;  called  normal  and  final  mode.  In  normal  mode 
one  is  using  the  code  to  predict  the  behavior  of  the  facility  based  on  ideal  test  conditions.  One 
provides  as  input  the  desired  operating  conditions:  final  temperature,  pressure  and  specific  heat 
ratio,  and  the  program  calculates  design  conditions  (based  on  the  test  turbine),  gas  properties  etc. 
This  mode  of  operation  can  exist  completely  outside  the  DAS  context  (i.e.  it  does  not  require  any 
DAS  format  file  input)  and  is  used  to  develop  initial  testing  conditions.  The  second  mode  of 
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operation  is  called  the  final  mode  and  it  differs  from  the  normal  mode  of  operation  in  that  it  takes 
gas  properties  recorded  during  the  fill  process  and  standards  reduction  process  and  generates  the 
test  gas  properties. 

Determining  the  Test  Gas  Properties 

In  normal  mode,  the  facility  model  assumes  a  desired  test  specific  heat  ratio  (y),  final  test 
temperature,  final  test  pressure,  and  the  two  types  of  gas  used  in  the  test  (x  for  the  first  component 
and  y  for  the  second  one).  The  program  calculates  mass  fractions  based  on  the  following  formula; 


M  .  _  Cy.x  (y-Tx) 

iVl-  — /  ^ 


mv 


'v,y 


(Yv-Y) 


(1) 


where  y  is  the  desired  test  condition  and  all  other  properties  are  functions  of  the  individual  gases  at 
the  present  temperature^.  All  properties  of  the  test  gas  mixture  are  functions  of  the  mass  ratio  and 
the  property  ratios  (gas  constant,  specific  heat,  etc.)  of  the  two  gases.  In  the  final  configuration 
mode,  one  takes  the  fill  data  (which  is  partial  pressure  and  temperature  data  for  the  initial  gas  and 
the  final  conditions),  and  uses  that  to  produce  a  mass  ratio. 
m„  .  . . _ 
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where  the  subscript  r  denotes  a  ratio  of  the  property  of  gas  y  to  gas  x.  Once  the  actual  mass  ratio 
has  been  obtained,  the  gas  property  calculation  proceeds  as  listed  in  the  normal  case.  One  can  see 
that  in  normal  mode,  one  finds  the  desired  test  conditions;  final  mode  takes  the  data  from  the  ran 
and  generates  the  test  gas  conditions. 

The  partial  pressures  and  temperatures  are  obtained  in  the  process  FCP.  This  process 
marks  the  partial  pressures  and  temperatures,  as  well  as  calibrates  sensors,  and  stores  these  values 
in  the  RCDF.  There  is  an  assumption  that  the  test  gas  properties  are  determined  just  from  these 
ratios.  However,  there  may  be  some  difference  between  the  final  fill  point  conditions  and  those  at 
which  the  facility  is  ran.  This  could  arise  due  to  changes  in  temperature,  venting  of  the  gas,  or  just 
leakage  over  time.  Thus  one  more  piece  of  information  is  needed  before  the  facility  model  can  be 
ran  in  final  mode,  that  is  the  initial  pressures  and  temperatures  in  the  supply  tank  and  the  test 
section  just  before  a  test. 

This  information  is  obtained  by  the  DAS  system  in  the  time  after  the  main  valve  is 
triggered,  but  before  the  main  valve  opens.  The  initial  conditions  are  determined  using  the 
"standards"  process.  In  this  process,  the  standards  channels  are  used  to  create  four  files  (time 


^  The  present  facility  model  only  models  gas  properties  as  a  function  of  temperature.  The  Calspan  analysis  shows 
how  this  data  could  be  modified  to  account  for  pressure  variations  as  well.  All  of  this  analysis  is  discussed  in  more 
detail  in  that  document 
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histories),  and  four  scalar  values.  The  four  files  are  the  defined  pressure  and  temperature  of  the 
supply  tank  and  test  section.  The  four  scalar  values  are  the  average  of  these  files  over  the 
beginning  of  the  test.  The  four  scalar  values,  and  the  time  range  over  which  they  were  averaged 
are  stored  in  the  RCDF  by  the  standards  process.  The  standards  process  controls  which 
instruments  are  used  to  create  these  files. 


Once  standards  has  been  run,  the  RCDF  is  sufficiently  complete  for  facility  to  run  in  final 
mode  and  generate  the  actual  gas  properties.  Thus  the  programs  which  must  be  run  to  determine 
the  gas  properties  used  in  the  test  are; 


Used  to  check  operation  of 
facility  and  determine 
loading  conditions 


Used  to  calibrate  sensors 
before  and  after  a  run  and 
to  mark  the  partial  pressures 
and  temperatures  at  the 
two  fill  points 

Used  to  acquire  data 
and  generate  DAS  files 


Used  to  reduce  standard 
channels  and  generate 
intial  conditions  in  supply 
tank  and  test-section 

Uses  fill  conditions 
to  determine  test  gas 
properties 


Operation  of  "mm" 

The  operation  of  the  program  "facility"  is  relatively  straight  forward.  However,  the 
operation  of  its'  X-1 1  interface  may  not  be  in  the  beginning  as  clear.  When  one  invokes  mm  there 
are  two  lines  of  options  (see  figure  1).  The  two  top  left  buttons  (full  printout,  and  startup 
dynamics)  control  the  number  of  data  points  plotted  in  the  output  files.  Basically,  if  one  would  like 
to  examine  the  main  valve  motion,  these  should  be  highlighted.  These  responses  are  put  into  the 
parameters  file.  The  next  three  buttons  control  file  operations,  running  the  model,  and  plotting  the 
results.  The  next  button  is  an  on-line  help  menu  and  the  last  button  terminates  the  program.  The 
buttons  on  the  second  layer  control  the  parameter  file  and  configuration  file  (except  for  the  button 
marked  "Final  Mode  Params".  The  far  left  button  (Parameters)  displays  all  the  values  which  go 
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into  the  parameters  files  (these  are  the  values  displayed  in  figure  1).  The  last  three  buttons  display 
all  the  values  that  go  into  the  configuration  file  (figure  2, 3,  and  4).  These  can  be  viewed  and 
edited. 

If  one  were  to  run  in  initial  mode,  one  could  go  to  the  file  options  button  and  load  in  an  old 
parameters  file,  or  one  could  use  the  X-1 1  information  to  create  a  parameters  file.  The  same  holds 
tme  for  the  configuration  file.  The  last  file  which  needs  to  be  specified  is  the  valve  data  file,  which 
could  either  be  an  old  or  a  new  file.  If  one  makes  any  changes  to  either  the  parameter  file  or  the 
configuration  file  after  they  have  been  loaded  in,  the  program  automatically  stores  these  changes  to 
temporary  files,  which  are  used  to  run  the  program  facility.  However,  if  one  wants  to  save  these 
files  for  future  use,  one  has  to  explicitly  save  them  using  the  file  options  button.  The  other  two 
options,  load  and  save  RCDF  are  used  only  if  the  program  is  being  run  in  final  configuration 
mode. 

If  one  wants  to  run  mm  in  final  configuration  mode,  one  would  load  or  create  parameter 
and  configuration  files  as  noted  above,  and  select  a  valve  data  file.  Then  one  could  go  into  the  final 
mode  set-up  by  pushing  the  "Final  Mode  Params"  button.  This  creates  a  window  (figure  5)  that 
has  all  of  the  parameters  set  by  "facility"  when  run  in  the  final  mode.  If  at  this  point  one  loads  a 
RCDF  (which  has  come  from  a  run,  and  has  been  operated  on  by  "standards")  one  should  see  that 
only  the  first  four  parameters  (the  partial  pressure  data)  have  values,  as  well  as  the  four  initial 
conditions.  One  operates  in  final  mode  by  pushing  the  "set  final  mode  button"  and  then  using  the 
run  model  button.  After  the  program  has  stopped,  the  values  for  the  test  gas  properties  will  be 
stored  into  a  temporary  RCDF  and  can  be  viewed,  but  not  edited,  in  the  final  mode  window.  If 
these  values  are  correct,  which  they  should  be,  one  then  pushes  the  "save  RCDF  file"  button  from 
the  file  options  menu.  This  will  result  in  a  statement  saying  that  if  one  continues,  these  values  will 
be  stored  in  the  RCDF  and  resulting  changes  in  the  RDRs  will  be  translated  to  all  data  files  in  the 
run.  Once  this  has  been  done,  the  action  is  also  logged  in  the  run  log  file. 

The  results  of  the  model  run  can  be  displayed  by  running  the  plot  mode  while  in  mm, 
running  "bdplot"  outside  of  mm,  or  by  using  "drp"  to  plot  the  DAS  files.  It  is  important  to  note 
that  while  one  can  be  in  any  directory  to  run  the  initial  mode  configuration,  one  should  run  the  final 
mode  configuration  from  the  run  directory  since  the  program  requires  an  RCDF,  broadcasts  the 
changes  in  the  RDR  to  all  existing  DAS  files  in  the  run  directory  and  either  creates  or  updates  the 
run  log  file. 

Future  Work: 

There  is  some  possibilities  of  improving  the  facility  model.  One  item  would  be  to  have  the 
facility  model  generate  the  partial  pressure  fill  conditions  when  running  in  initial  mode.  In 
addition,  the  gas  properties  presently  used  in  the  model  are  ideal  gas  properties  and  are  only  a 
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function  of  temperature.  It  may  be  desirable  to  updated  the  model  to  include  pressure  dependency 
which  would  require  some  rework  of  the  gas  selection  subroutines  and  should  probably  be 
evaluated  in  the  context  of  an  overall  uncertainty  analysis  of  the  facility.  Clearly  one  of  the  main 
goals  of  running  the  facility  model  in  normal  mode  is  deciding  how  to  set  the  opening  and  closing 
reservoirs  on  the  main  valve.  The  issues  involved  with  the  six  chamber  model  need  to  be  resolved. 
Once  that  is  done,  the  facility  program  will  have  to  be  updated  to  reflect  these  changes. 
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Current  Data  Group: 


Paraiiieters 


Sprin:^i^d  iS^e 


Figure  1:  X-11  Interface  Screen  for  "facility"  (mm] 
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Figure  2:  Volumes,  Orifice,  Valve  Windovv' 
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>la&aibdel 


Full  PHutoit 


File 


Riiti  Facility  HodalifFlot 


Current  Data  Croup: 


Par I  IVol  umes  ^  Or i  f i  ce ,  Val 


Actuator  Cylinder  and  Gas|| 


Sp  rl  n9;;;ari^s  Brake 


I 


y,NV^^%S%XV^\VVV»>ywN\VK\Vw\VS\\VW%\VV\^WSSVy^VW^WV^SVVNVWWWWS\'^y^VVWV^^^NywW^ 


Bore  Dianeter  (inches)  |p>0 


Rod  Dianeter  (inches)  l0*0 


Nunber  of  Rods 


Mflx  Actuator  Stroke  (inches)  [p*0 


Piston  Thickness  (inches)  [p«0 _ 

>  Stroke  Offset,  S  Q  X=0  (inches)  |p,0  \m;i 

I  Offset  to  Start  of  Transfer  Port  (inches)  [p>0 
'■■  End  of  Transfer  Port  (inches)  [p.O 
Transfer  Port  Area  (sq  in)  [p,0  | 


;  Actuator  Fixed  Bleed  Area  (sq  in)  [p^O  I? 


Leakage  Around  Pison  (sq  in)  [p>0 _ 


Suitch  Position  (inches)  LO.O 


\N%NS'^VHSNSVW'W\S\\VS\\KSSN\\\\\V\S'N*N\S\\\\NVwSS\\NNS\V\%SNX\\\\\VV\\NS\S**V\SS\NN%\SVN^N\SSS\\\V\\SN\N\\NVVWVNV\NSNN\VS.VWW 


Opening  Gas  Supply  Volume  (cubic  H)  ip.Q 


I  Opening  Line  Volume  (cubic  in)  [0*0 


Closing  Line  Volune  (cubic  in)  10*0 _ 


Closing  Cas  Supply  Volume  (cubic  H)  ip*o 

iViYiYiViYuiYiYiYiViYiYM^^^^^  . 

Vent  Indicator  mure»j[' 


:  Opening  Cas  Supply  Orifice  (sq  in)  [p*0  i 

;  Closing  Cas  Supply  Orifice  (sq  in)  [p.O 
■  Cylinder  Opening  End  Vent  Orifice  (sq  in)  [ipTo 
i  Cylinder  Closing  End  Vent  Orifice  (sq  in)  [pVo 
:  Dead  Volume  in  Cylinder  Head  (cubic  in)  (p.O 


,SNV\NN\\S\\SVS\N\\\VIiNN\S\\\S\NVVS\VN%S\NV>\V\SS\!V^’ 


<f  ''  '  J‘  f  r 

Figure  3:  Actuator  Cylinder  and  Gas  Window 
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Figure  4:  Spring  and  Brake 
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Current  Data  Group: 


final  Mode  Parans 


Valunes.  Qpii'^ce,  Va^ii^[tetHator  Cylinder  ati<t  Gas i:  Spring  and  Srakei 


Partial  Pressure  of  Gas  #1  [kPa]  6 

loo  ^ 

Tenperature  of  Gas  il  [k]  oToo 

— 

ji:  Final  Fill  Pressure  [kPa]  0.00 

5  :•:  _ _ 

Final  Fill  Tenperature  [kl  0.00 

- 

5:;  Test  Gas  hole  Frac  1  [p.OO 

5  ;  -wxX.*  vx^  »  •W/'  .'xVs*'"' V*  j'  ■>  ^  j-  x 

||i  Test  Gas  Mole  Frac  2  [p.OO  |s® 

Molecular  Ueight  [p.OO  i 

11:  Specific  Heat  Coef  1  [p.OO 

j  ;.v.%v.%v.v.%vA«-v.%%WwvAWA'^V4-Ar.W;^A^^^V4Vw^AW*VAV.w«W:/A»Af.*.wj*.%w.*>-* 

*; 

li:  Specific  Heat  Coef  2  [p.OO 

;  ;.V,%V.%V.W.W.WA-A%W,VAV-WAWJWfc/AWA"A%*<WW^VA*«WAr--rfVV‘.P-*A-.%%-AV.WA*:*:* 

Specific  Heat  Coef  3  [p.OO  \m 

\  ;/.\-A%%%%%-AWA^---WAfl.WAVWiaArAAAVAVW^VWA*i^A'rfVVWWA\WlAW.Vrf'^AV.V.-.-.VAV;*:- 

ii:  Spec  Heat  Ratio  1  [p.OO 

y  . .  .  :  i 

!**  if  *  'J; 

•¥ 

i|:  Spec  Heat  Ratio  2  [p.OO  \-- 

i|;  Spec  Heat  Ratio  3  [p.OO  |  ^ 

!  V  . 

Ii  Prandtl  Nunber  1  [p.OO  |  - 

*  •  *  •  .<• 

Prandtl  Nunber  2  [p.OO  |®: 

!  Prandtl  Nunber  3  [p.OO  I 

tIjRal  P4o^  Parswwt»i*s  <coRt:»)i 


^  Bynanic  Viscosity  1  [p,00 


I  Dynanic  Viscosity  2  10*00 


I  Dynanic  Viscosity  3  lO.OO 


I  Gas  Constant  (Ru/MU)  [p.OO 


Coadltioi^s 

. "»>W-jywy*V»V0S.» . <»yy.yrw»n»w 


5  Upstrean  Des  Stag  Tenp  (deg  K)  [p*00 


Upstrean  Des  Stag  Press  (kPa)  [0*00 


Total  Des.  Ueight  Flou  (kg/sec)  [p*00 


Des.  Ueight  Flou  thru  Turb  <kg/sec>  [p,0Q 


Design  Prandtl  Nunber  of  Gas  [p.OO 


Ssata 


Supply  Tank  Tenp  (deg  K)  [p.OO 


Supply  Tank  Press  <kPa>  [p.OO 


Test  Section  Tenp  (deg  K)  [O.QQ 


>  Test  Section  Press  (kPa)  [p.OO 
[$  Standards  Int  Start  [O.QQ 


Standards  Int  End  [p  .QQ 


. . . . 


■i, 

::l^ 


Figure  5:  Final  Mode  Parameters  Window 


Original  Math  Model  Description 


by:  Mike  Stanton 
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1  Gas  Properties 


The  thermodynamic  properties  of  a  gas  mixture  are  calculated  from  the 
properties  of  the  component  gases.  A  one  or  two  component  gas  can  be  se¬ 
lected.  The  gases  in  the  database  are:  Air.  Argon,  Freoni2,  N2,  and  CO^- 
The  gas  component  concentrations  are  those  whose  mixture  specific  heat  ra¬ 
tio  give  the  best  match  to  the  desired  specific  heat  ratio.  The  database  for 
the  the  gas  components  was  derived  from  the  following  sources: 

1.  Hilsenrath,  NBS  CIR564,  QC286.H655  -  all  but  Freonx^- 

2.  Ashrae,  Thermophysical  Properties  of  Refrigerants,  TP492.A54  -  for 
Freon\2‘ 


1.1  List  of  Variables 


Cp.TTlXr 

specific  heat  at  constant  pressure  of  the  gas  mixture 

J 

kg^^K 

specific  heat  at  constant  pressure  of  gas  component  i 

J 

kg -  OK 

specific  heat  at  constant  volume  of  the  gas  mixture 

J 

^vA 

specific  heat  at  constant  volume  of  gas  component  i 

J 

kg^oK 

Tnitx 

ratio  of  specific  heats  of  the  gas  mixture 

li 

ratio  of  specific  heats  of  gas  component  i 

\V 

kmix 

thermal  conductivity  of  the  gas  mixture 

k: 

thermal  conductivity  of  gas  component  i 

\v 

m-°K 

l^mix 

dynamic  viscosity  of  the  gas  mixture 

dynamic  viscosity  of  gas  component  i 

*77?” 

Frjnix 

Prandtl  Number  of  the  gas  mixture 

moleculax  weight  of  the  gas  mixture 

moie 

Mi 

Ro 

molecular  weight  of  gas  component  i 

Gas  Constant 

Q 

m^lc 

moic-®A“ 

Rgaa 

gas  constant  for  the  mixture 

J 

kg-oK 

Xi 

mole  fraction  of  gas  component  i 

Xi 

concentration  of  gas  component  i 

qrams 

gram  of  mix 

13 
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Fdcility  Configuration  Model 


1.2  Defining  Relationships  for  Properties  of  the  Gas  Mixture 
The  ga5  mixture  properties  are  calculated  via  the  following: 

^  ^  _ 72Cp.i(7-7l) _  (1) 

^  72Cp,i(7-7:)  -  7iCp,2(7-72) 

subject  to  X2  =  Miiumuin(l.wY2)  and  X2  =  MaximumfO.jXj) 

Xi  =  1.0 -.Y2  (2) 


(3) 

(4) 

(5) 

(6) 

(7) 

(8) 
(9) 

(10) 

(11) 

(12) 

(13) 

(14) 
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Facility  Configuration  Model 


2  Facility  Geometric  Parameters 

These  consist  of  volumes,  flow  areas,  and  pressure  face  areas. 


Volumes: 


Vb 

blowdown  supply  tank 

fixed 

turbine  inlet  passage 

fixed 

14 

turbine  exhaust  passage 

fixed 

14 

down  stream  of  throttle 

fixed 

V, 

dump  tank 

fixed 

Vs 

actuator  reservoir  (opening) 

fixed 

Vs 

actuator  reservoir  (closing) 

fixed 

Vs 

actuator  bore  (opening) 

variable 

V7 

actuator  bore  (closing) 

variable 

Vf6 

actuator  line  (opening) 

fcced 

Vf7 

actuator  line  (closing) 

fixed 

Flow  Areas: 

.4oi 

-4i2 

•4.14 

-423 

-434 

-456 

.487 

Atv 

Aev 

Af67 

Av67 

-467 


main  valve 

variable 

turbine  throat 

fixed 

boundary  layer  bleed  =  0.34i2 

fixed 

iTTn  P'y 

throttle  area  =  AiiJ-rj^j^ 

variable 

dump  tank  inlet 

fixed 

actuator  supply  (opening) 

fixed 

actuator  supply  (closing) 

fixed 

actuator  vent  (opening) 

fixed 

actuator  vent  (closing) 

fixed 

actuator  fixed  leakage 

fixed 

actuator  transfer  port 

variable 

actuator  leakage  =  4f67  +  4^67 

variable 

Pressure  Face  Areas: 

Ao  area  of  actuator  piston  (opening  face)  fixed 
4cr  area  of  actuator  piston  (closing  face)  fixed 


Facility  CoaUguration  Model 


3  Facility  Chamber  Pressures  and  Temperatures 

Pressures: 

po  blowdown  supply  tank 
Pi  turbine  inlet  passage 
P2  turbine  exhaust  passage 
P2  down  stream  of  throttle 
P4  dump  tank 

ps  actuator  reservoir  (opening) 
ps  actuator  reservoir  (closing) 
pe  actuator  bore  (opening) 
pr  actuator  bore  (closing) 

Temperatures: 

To  blowdown  supply  tank 
Ti  turbine  inlet  passage 
T2  turbine  exhaust  passage 

T3  down  stream  of  throttle 

r4  dump  tank 

Ts  actuator  reservoir  (opening) 

Ts  actuator  reservoir  (closing) 

Te  actuator  bore  (opening) 

Tt  actuator  bore  (closing) 

r,j  temperature  of  the  upstream  gas 
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4  Flow  Rates  at  Orifices 

The  flow  rates  are  calculated  for  adiabatic  flows.  The  flow  conditions 
are  considered  to  be  either  choked  or  unchoked  flow. 


4.1  List  of  Variables 


Pi 

Pi 

Aij 

Wij 

w 


P^p 

Pdown 


P  critical 

'^doxun 

Tii 


pressure,  chamber  i 

pressure  in  adjacent  chamber,  chamber  j 
area  of  the  orifice 

flow  rate  through  the  orince  from  chamber  i  to  chamber  j 
absolute  value  of  flow  rate 

chamber  pressure,  high  pressure  side  =  Maximum(p:,pj) 

chamber  pressure,  low  pressure  side  =  Minimuin(p{,pj) 

critical  pressure  at  which  the  flow  is  choked 

chamber  temperature,  high  pressure  side 

chamber  temperature,  low  pressure  side 

chamber  temperature,  high  pressure  side 

ratio  of  specific  heats  of  the  gas  mixture  =  7m:*r 


4.2  Defining  Relationships  for  Orifice  Flow  Rates 

/  2  - 
P  critical  “  ^  ^  ^ 

a  Pdorvn  <  Pcritia^u  flow  is  clioked  and: 


_  /  /  ^  Pup-^ij 

""  ~  V^^7  +  l  ■ 

else  if  Pdown  >  PcTiticah  tbe  flow  is  uncboked  and: 


w 


(15) 


(16) 


J(Jl-)(l-(g^^)-(^)^  (17) 

y  7  —  1  Pup  Pup  y  ^gasd-up 


W{j  = 


Pi-Pj 

\pi-pi\ 


'It; 


(18) 
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Facility  Configurar/on  Model 


5  Turbine  and  Brake  Properties 

5.1  List  of  Variables 


Pi 

P2 

Pratio 

TTi 

TT2 

T'ratio 

Mz 

TTi 

T 

■‘■room 
'^gaa  to  wall 

ht 

Tt 

% 

Ir 

W5 

Re 

Dtip 

% 

%r 

hb 

Qh 

h 

^b,max 

'Rdrwm 

“^dTum 

^p,dTiim 

B 

Br 

Be 

Bo 

Bi 

B2 


t 

^brake  on 
^brake  off 


turbine  inlet  pressure 
turbine  exit  pressure 
turbine  pressure  ratio 
turbine  inlet  temperature 
turbine  exit  temperature 
turbine  temperature  ratio 
boundary  layer  bleed 

time  rate  of  change  of  turbine  inlet  temperature 

ambient  temperature 

turbine  gas  to  wall  temperature 

adiabatic  efnciency  of  the  turbine 

turbine  power 

turbine  shaft  torque 

net  shaft  torque 

rotor  moment  of  inertia 

turbine  shaft  speed 

scaled  turbine  shaft  speed 

turbine  shaft  design  speed 

scaled  turbine  shaft  design  speed 

turbine  Reynolds  Number 

tip  diameter  for  turbine  rotor 

actual  brake  torque 

required  brake  torque 

brake  power 

total  energy  absorbed  by  the  brake 
brake  coefficient 

shaft  speed  at  which  maximum  brake  torque  develops 

brake  temperature  (with  no  dissipation) 

mass  of  brake  drum 

heat  capacity  of  brake  drum 

actual  brake  excitation 

required  brake  excitation 

total  brake  constant 

first  order  brake  constant 

second  order  brake  constant 

third  order  brake  constant 

brake  time  constant 

run  time 

time  at  which  braking  action  starts 
time  at  which  braking  action  stops 
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5.2  Defining  Relationships  for  Turbine  and  Brake  Properties 

P2 

Praiio  _ 

Pi 

(20) 

Tratio  =  1  —  €a(l  —  Pratio  ) 

(21) 

TTi  =  Ti 

(22) 

TT2  =  TTi  -Tratio 

(23) 

y/Traiio 
—  A\2  •  ■"  p 

* ratio 

(24) 

ht  =  Wi2-Cp,mix'{TTi  -TT2) 

(25) 

II 

(26) 

u;. 

(27) 

^'ds 

<^sds  -  ^ 

(28) 

Re  = 

-UtipP^mtx 

(29) 

' 

TTi 

'^gas  to  wall  “  nr» 

-iroom 

(30) 

TTlIr'*^sds 

~  zz~  '2,/m 

(31) 

^l  kl^is 

(32) 

if  t  <  tbrake  on  OT  t  >  thrakt  off 

0 

II 

(33) 

else 

if  excitation  is  not  preprogTaTmned  then 

B  =  Br 

(34) 
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I 


otherwise  for  a  preprogrammed  excitation 

Be  —  Bq  +  Bi{t  —  tiyrake  on)  “T  B2{t  —  hrake.  on) 


and 


t—thrnk^  on  V 

B  =  5c-(l-e - r) 


hi  =  ki-B^ 


W5 


^  '  ''^b,Tnax^ 


Wj 

T„  =  %-% 

rp  _  _ Qb  ,  rp 

^  drum  ~  ^  ^  room 

T^drum  '  ^p.drum 


(35) 

(36) 

(37) 

(38) 

(39) 

(40) 


6  Main  Slider  Valve  and  Actuator  Variables 


6.1  List  of  Variables 


I 

V 

Ts 

^ksc 


^kso 

So 

9 

Ft 

F,t  . 
Fp 
Fs 
Fa 

Ft 

K,c 

Kso 


sliding 

l^atatic 

m 


9 


the  disDlacemeat  of  main  slider  valve  at  which  the  front  seal  is  fully  compressed 

velocity  of  main  slider  valve  =  velocity  of  actuator  piston 

radius  from  centerline  of  chamber  1  to  the  back  seal 

radius  from  centerline  of  chamber  1  to  the  front  seal 

displacement  at  which  the  front  seal  is  closed  but  uncompressed 

X  displacement  at  which  the  front  slider  spring  is  in  equilibrium 

X  displacement  at  which  the  back  slider  spring  is  in  equilibrium 

offset  of  actuator  piston  from  front  of  bore  at  x=0 

maximum  area  of  main  valve 

angle  from  vertical  of  the  main  valve  bevel  • 

total  force  on  slider  assembly 

total  force  excepting  friction  on  slider  assembly 

pressure  force  contribution  to  total  force  on  slider  assembly 

spring  force  contribution  to  total  force  on  slider  assembly 

actuator  force  contribution  to  total  force  on  slider  assembly 

frictional  force  contribution  to  total  force  on  slider  assembly 

spring  constant  for  front  spring 

spring  constant  for  back  spring 

coefficient  of  sliding  friction  for  the  main  slider  valve 

coefficient  of  static  friction  for  the  main  slider  valve 

mass  of  sliding  assembly 

gravitational  acceleration 


I 

I 
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6.2  Defining  Relationships  for  Main  Valve  and  Actuator 


^01  —  0  X  ^  Xjji 

.4oi  =  A  =  2-rm(s-Xjn)sind-r-(x—Xm)^sm^0cosd  x  >  x-nO-ndA  <  .4m 

(42) 

i4oi  =  .4m  2  >  2m  and  A  >  Am  (43) 

=  .4o(5o  +  2) -f  Vps  (44) 

V-r  —  Ac^So  +  2)  -r  VfT  ■  ;  (45) 

Fp  =  -T(rm^ -r,^)(po -Pi)  X  <  Xm  (46) 

Fp  =  0  2  >  im  (47) 

Fa  =  P6  •  -4o  -  P~  •  -4c  (48) 

F,  =  -Ksc{x  -  Xksc)  X  <  Xksc  (49) 

Fj  —  0  Xksc  —  X  Xkso  (80) 

Fj  —  —-^50(2  —  Xkso)  X  ^  Xkso  (^4) 

Fst  =  Fp  +  Fa-rFs  (52) 

Ff  =  PslidmgTTig  V  #  0  (53) 

~  ~  I  ■  Pstaticnag  V  =  0.  [Fsil  ">  fistatic^ag  (54) 

\Xst\ 

Fj  =  -Fst  V  =  Q,  IFstl  <  fistaticTng  (55) 

Ft  =  Fst  +  Ff  (56) 


7  Method  of  Integration 

7.1  List  of  State  Variables 

Us  shaft  speed 

V  slider  velocity 

X  slider  displacement 

mi  mass  of  gas  in  chamber  i  (i=0...8) 

Ei  internal  energy  of  gas  in  chamber  i  (i=0...8) 

TT\  turbine  inlet  gas  temperature 
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7.2  Calculation  of  the  Derivatives  of  State  Variables 

Tn 

m 

V 


V  = 

X  = 

TTIq  =  -lOQi 

mi  =  +  Woi  -  “^^12  -  '^14 
1712  =  “r^l2  -  ^^23 
m3  =  -f  1^23  - 

7714  =  -ry^u  “T  ^^34  "T  y^sv  +  y^TV 

m^  —  ^^*56 

ms  =  -r  iL'se  -  y^S7  -  y^ev 
ihr  =  -ru;67 -r 

Thg  =  —  71787 

^0  ~  ^01  ^p,m:x^01 

£1  =  +7^01^p,mxz^l  “  '^12^p^mtr^l2  'I^14^p,mtx^l4 

J&2  “  "r^l2^p,Tn:x^l2  ^23^,mix^23 

£3  =  +W2ZCp^mixT2Z  -  yyZ4Cp,mixTz4 
?4  “  +‘U;i4Cp,;ntxri4  +  'li^34Cp.7n:x?34  +  V^'eVCp^mixTsV  +  'iVTVCp^mixTTV 

£5  =  -W56Cp,^n6 

JSg  =  +t/;5eCp^^^T56  —  y^erCp^TnicTer  ^  y^evCp.mix'I'ev 

(uy 

Er  =  +We7Cp^mixTsr  +  yyS7Cp,mixTs7  -  WrVCp^mixTTV 

Uit- 

Es  ^  — y^S7^p,7^i^^S7 
UiS 

El  —  CvTuixTiihi 

TT,  =  ^  :—— 

he  integration  scheme  is 
ri+M 


772i  Cy^ 


(57) 

(58) 

(59) 

(60) 
(61) 
(62) 

(63) 

(64) 

(65) 

(66) 

(67) 

(68) 

(69) 

(70) 
(Tl) 

(72) 

(73) 

(74) 

(75) 

(76) 
(T7) 

(78) 


,mxx 


rt+At 

dS  =  St+M-St  =  St  At 


(79) 
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•^t+At  —  St -r  St^t  (SO) 

Pressure  and  T  emperature  are  updated  from  the  integrated  state  variables  : 

Ti  = 

Pi  = 

8  User  Specified  Parameters  to  the  Facility 

Model  Program 

IPNTOUT  1  =  full  printout 
0  =  short  printout 

ISD  1  =  startup  dynamics  are  modelled 

0  =  start  with  turbine  inlet  and  outlet  at  design  pressure  and  temper¬ 
ature  ratio.  Also  inlet  conditions  =  supply  conditions 

ITURB  0  =  user  specified  turbine 
3  =  T3  Turbine 

7  =  MIT  Model  Ace  Turbine 

8  =  Full  Scale  Ace  Turbine 

PIDP  (entered  when  ITURB=0)  design  outlet  pressure  to  inlet  pressure 
ratio 

ETADP  (entered  when  ITURB=0)  design  adiabatic  efficiency 
GDP  (entered  when  ITURB=0)  design  7 
REDP  (entered  when  ITURB =0)  design  Reynold’s  number 
AT  (entered  when  ITURB =0)  turbine  throat  area  (in^) 

DTIP  (entered  when  ITURB=0)  rotor  diameter  (in) 

TGTWDP  (entered  when  ITURB=0)  design  gas  to  wall  temperature 
IR  (entered  when  ITURB=0)  rotor  moment  of  inertia  {kg  —  m?) 

NSO  design  operating  speed  (RPM) 


Ei 

(81) 

EiRgas 

(82) 

V{Cy^gcis 
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IBEX  1  =  brake  excitation  is  computed 

2  =  brake  excitation  is  preprogrammed  via  tbe  quadratic  constant,  Be, 
Be  =  Bq  +  B\t  T  B-^t^ 

BO  (entered  when  IBEX=2)  zero  order  constant 
Bl  (entered  when  IBEX=2)  first  order  constant 
B2  (entered  when  IBEX=2)  second  order  constant 
POPSI  initial  supply  tank  pressure  (psia) 

TOIN  initial  supply  tank  temperature  where  a  positive  number  indicates  a 
temperature  in  °K  and  a  negative  number  in  ‘’P 

P4PSI  initial  dump  tank  pressure  (psia) 

T4IN  initial  dump  tank  temperature  where  a  positive  number  indicates  a 
temperature  in  °K  and  a  negative  number  in  °F 

GD  desired  7  for  the  test  gas 

NCMP  number  of  gas  components  in  the  gas  mixture  (  1  or  2  ) 

ICMP  index  of  gas  components  in  the  gas  mixture 

1  for  .4ir 

2  for  Argon 

3  for  Freoni2 

4  for  N2 

5  for  CO2 

P5PSI  actuator  reservoir  (opening)  initial  pressure  (psia) 

P8PSI  actuator  reservoir  (closing)  initial  pressure  (psia) 

PVR.PSI  vent  reservoir  initial  pressure  (psia) 

PVAOPSI  actuator  bore  initial  pressure  (psia) 

TSIM  length  of  simulation  (seconds) 

NPMS  number  of  time  steps  per  millisecond  -  10  is  recommended.  MIT 
Experience  has  shown  that  too  few  time  steps  per  millisecond  will  show 
up  as  pressure  oscillations  on  the  order  of  a  time  step  in  the  solution. 
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TOPEN  actuator  valve  opening  time  (seconds) 

TCLOSE  actuator  valve  closing  time  (seconds) 

TON  brake  turn  on  time  (seconds) 

TOFF  brake  turn  off  time  (seconds) 

9  Facility  Configuration  File  Parameters 

VO  SUPPLY  TANK 
VI  TURBINE  INLET 
V2  TURBINE  EXHAUST 
V3  DSV  INLET 
V4  DUMP  TANK 

ADSO  EXIT  PIPES  TO  DUMP  (in-) 

ISS  =1  FOR  SLIDING  SLEEVE  (ATAR),  =0  FOR  COWCATCHER  (MT) 

RO  VALVE  ANNULUS,  OUTER  RADIUS  (INCHES) 

RI  VALVE  ANNULUS,  INNER  RADIUS  (INCHES) 

RM  MAIN  SEAL  RADIUS  (INCHES) 

XM  MAIN  SEAL  BREAKAWAY  (INCHES) 

RS  BACK  SEAL  BREAKAWAY  (INCHES) 

XS  BACK  SEAL  BREAKAWAY  (INCHES) 

THETA  ENTRANCE  ANGLE  FROM  VERTICAL  (DEGREES) 

MASS  TOTAL  MOVING  MASS  (KG) 

XKSC  NO  FORCE  POSITION,  MAIN  VALVE  CLOSING  SPRING  (INCHES) 

KSC  SPRING  CONSTANT  (7^),  CLOSING  SPRING 

XKSO  NO  FORCE  POSITION,  VALVE  OPENING  SPRING  (INCHES) 
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KSO  SPRING  CONSTANT  (7^),  CLOSING  SPRING 
CSF  COEFFICIENT  OF  STATIC  FRICTION 
CDF  COEFFICIENT  OF  SLIDING  FRICTION 
DB  BORE  DLA.METER  (mCHES) 

DR  ROD  DLA.METER  (INCHES) 

NRODS  2  FOR  ATAR,  1  FOR  xMIT 

SMAX  MAXIMUM  ACTUATOR  STROKE  (INCHES) 

SP  PISTON  THICKNESS  (INCHES) 

50  STROKE  OFFSET,  S  @  X=0  (INCHES) 

51  OFFSET  TO  START  OF  TRANSFER  PORT  (INCHES) 

52  END  OF  TRANSFER  PORT  (INCHES) 

AV67  TRANSFER  PORT  AREA  {in^) 

AF67  ACTUATOR  FIXED  BLEED  AREA  {in‘) 

Vo  OPENING  GAS  SUPPLY  VOLUME  (m^) 

VF6  OPENING  LINE  VOLUME  (m^) 

VF7  CLOSING  LINE  VOLUME  (m^) 

V8  CLOSING  GAS  SUPPLY  VOLUME  (m^) 

IVNT  VENT  TO  DUMP=0,  VENT  TO  ROOM=1 
A056  OPENING  GAS  SUPPLY  ORIFICE  {in^) 

A07V  CYLINDER  CLOSING  END  VENT  ORIFICE  (tV) 
AC87  CLOSING  GAS  SUPPLY  ORIFICE  (in^) 

AC6V  CYLINDER  OPENING  END  VENT  ORIFICE  (in^) 
KB  BRAKE  CONSTANT 

NB  BRAKE  SPEED  @  MAXIMUM  TORQUE  (rpm) 
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MDRUM  BR.A.KE  DRUM  MASS  (KG) 

CPD  BBAKE  DRUM  SPECIFIC  HEAT  (7;^) 

TAUB  BRAKE  TIME  CONSTANT  (SECONDS) 

10  Operating  the  Facility  Model  Programs 

The  programs,  ”bd”  and  ”bdplt”,  are  the  programs  which  produce  the 
computational  results  of  the  modeling  and  the  plot.  The  program,  "bd”,  is 
operated  by  typing  ”bd”  and  answering  the  prompts  for  the  user  selected 
parameters.  The  facility  configuration  file  named  ” configuration.dat”  must 
exist  in  the  present  working  directory.  Two  output  files  are  produced  by 
this  program.  The  file,  "BD.OUT”  contains  a  listing  of  the  computational 
results  of  the  run.  The  file,  ’’BD.PLT”  contains  the  data  file  which  is  used 
as  input  to  "bdplt”. 

The  program,  ”bdplt”,  is  operated  by  typing  ”bdplt”.  This  program 
produces  the  plots.  Plots  can  be  generated  on  the  screen  or  on  a  laser 
plotter.  The  environment  variable,  mplotfile.  is  used  to  pass  the  names  of 
the  input  file  to  ”bdplt”.  The  variable  mplotfile  contains  the  name  of  the 
plotfile  produced  by  ”bd”.  It  is  set  by  typing  the  unix  command  "setenv 
mplotfile  BD.PLT”  or  whatever  else  the  input  file  is  named. 
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Program /flc//in’  has  two  basic  forms  of  supplying  input  data.  The  first  consists  of  a  file  contain¬ 
ing  facility  configuration  parameters.  There  are  seven  distinct  groups  of  data  which  are  entered 
using  the  FORTRAN  namelist  convention. 

NAMELIST  CONTENTS 
FDl  Facility  Volumes 

FD2  Flow  Path  Orifice 

VDl  Valve  Geometry  Variables 

VD2  Impact  Spring  Data 

VD3  Actuator  Cylinder  Geometry 

VD4  Actuator  Gas  Supply  Geometry 

BRK  Brake  Constants 

Each  namelist  has  one  or  more  parameter  names  which  are  used  to  enter  the  associated  input 
data.  The  namelist  parameters  are  presented  for  the  user’s  convenience: 

FDl  Facility  Volume  Data  (Cubic  Feet) 

VO  SUPPLY  TANTC 
VI  TURBINE  INXET 

V2  TURBINE  EXHAUST 

V3  DSV  INLET 
V4  DUMP  TANK 

FD2  Facility  Flow  Path  Orifice  (Square  Inches) 

ADS 0  EXIT  PIPES  TO  DUMP,  A34: 2-20"DIA,  CD=.62 

VDl  Valve  Geometry  Variables 

ISS  INTEGER  INDICATOR  [0=MIT  COWCATCHER;  1=ATARR  SLIDING  SLEEVE], 
RO  VALVE  ANNULUS,  OUTER  RADIUS  (INCHES) 

RI  VALVE  ANNULUS,  INNER  RADIUS  (INCHES) 

RM  MAIN  SEAL  RADIUS  (INCHES) 

XM  MAIN  SEAL  BREAKAWAY  (INCHES) 

RS  BACK  SEAL  RADIUS  (INCHES) 

XS  BACK  SEAL  BREAKAWAY  (INCHES) 

THETA  ENTRANCE  ANGLE  FROM  VERTICAL  (DEGREES) 

MASS  TOTAL  MOVING  MASS  (KG) 


V 
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VD2  Impact  Spring  Data 

XKSC  ZERO  FORCE  POSITION,  MAIN  VALVE  CLOSING  SPRING  (INCHES) 

KSC  SPRING  CONSTANT  [CLOSING  SPRING]  (LBF/IN) 

XKSO  ZERO  FORCE  POSITION  [VALVE  OPENING  SPRING]  (INCHES) 

KSO  SPRING  CONSTANT  [OPENING  SPRING]  (LBF/IN) 

CSF  COEFFICIENT  OF  STATIC  FRICTION 

COEFFICIENTS  OF  DYNAMIC  FRICTION  FOR  SLIDER  RANGES  SHOWN: 
CDl  OTOXl 
CD2  X1TOX2 
CD3  X2  TO  X3 

CD4  X3  TO  END 

SLIDER  POSITIONS  FOR  DYNAMIC  FRICTION  COEFRCIENTS: 

XI  LOCATION  1  (INCHES) 

X2  LOCATION  2  (INCHES) 

X3  LOCATION  3  (INCHES) 

VD3  Actuator  Cylinder  Geometry 

DB  BORE  DIAMETER  (INCHES) 

DR  ROD  DIAMETER  (INCHES) 

NRODS  NUMBER  OF  RODS  [  1  FOR  MIT;  2  FOR  ATARR] 

SMAX  MAXIMUM  ACTUATOR  STROKE  (INCHES) 

SP  PISTON  THICKNESS  (INCHES) 

50  STROKE  OFFSET,  S  @  X=0  (INCHES) 

5 1  OFFSET  TO  START  OF  TRANSFER  PORT  (INCHES) 

52  END  OF  TRANSFER  PORT  (INCHES) 

AV67  TRANSFER  PORT  AREA  (IN**2) 

AF67  ACTUATOR  FIXED  BLEED  AREA  aN**2) 

ALEAK  LEAKAGE  AROUND  PISTON  (IN**2) 

SSW  SWITCH  POSITION  (INCHES) 


^  - 
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VD4  Actuator  Gas  Supply  Geometry 

V5  OPENING  GAS  SUPPLY  VOLUME  (M^*3) 

VF6  OPENING  LINE  VOLUME  (IN**3) 

VF7  CLOSING  LINE  VOLUME  aN**3) 

V8  CLOSING  GAS  SUPPLY  VOLUME  (M**3) 

IVNT  VENT  INDICATOR  [0=DUMP;  l=ROOM] 

A056  OPENING  GAS  SUPPLY  ORIFICE  (IN*=^2) 

A07V  CLOSING  GAS  SUPPLY  ORIHCE  (IN*="2) 

AC87  CYLINDER  CLOSING  END  VENT  ORIFICE  (IN**2) 

AC6V  CYLINDER  OPENING  END  VENT  ORIHCE  (IN**2) 

VDEAD  DEAD  VOLUME  IN  CYLINDER  HEAD  (IN**3) 

BRK  Brake  Constants 

KB  BRAKE  CONSTANT  (W-S**2  /  TESLA  =^=^2) 

NB  BRAKE  SPEED  @  M.AXIMUM  TORQUE  (RPM) 

MDRUM  BRAKE  DRUM  MASS  (KG) 

CPD  BRAKE  DRUM  SPECIFIC  HEAT  (J/KG-K) 

TAUB  BRAKE  TIME  CONSTANT  (SECONDS) 

Guidelines  for  entering  each  namelist  name  and  parameters: 

First  line  contains  a  dollar  sign  ($)  in  column  2  and  the  namelist  name  starting  in  column  3.  For 
each  namelist  parameter,  open  a  new  line  and  enter  its  name  anywhere  beyond  column  1  fol¬ 
lowed  by  an  equal  sign  (=),  the  value  to  be  assigned  the  parameter,  and  a  terminating  comma  (,). 
For  example,  the  line  "  V3  =  250,"  appearing  in  the  FDl  namelist  input  stream  would  specify  a 
DSV  inlet  volume  of  250  cubic  feet.  To  terminate  input  for  a  namelist,  its  last  line  should  con¬ 
tain  the  four  characters  "SEND"  starting  in  column  2. 

The  second  form  of  input  to  program  facility  consists  of  the  responses  to  screen  prompts  for 
input  data  These  responses  may  be  entered  interactively  or  placed  in  an  ASCII  file  to  be  read  by 
the  program.  Since  the  data  are  order  dependent,  there  is  no  recovery  from  input  errors  when  a 
response  file  is  used.  Even  considering  this  limitation,  it  is  preferable  to  use  this  form  of  input 
because  it  provides  an  automatic  record  and  a  convenient  platform  for  changes.  For  the  most 
part,  a  response  file  consists  of  lines  of  single  numbers  whose  meanings  are  explained  in  the  out¬ 
line  presented  below. 
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Outline  of  Response  File  Entries  for  the  Facility  Model  Program 


Parameter 

IFRCDF 

IPNTOUT 

ISD 

ITURB 


NSO 

IBEX 


POPSI 

TOIN 

P4PSI 

T4IN 

GD 


Usage  in  Facility  Model  Program 

MODE  INDICATOR  [0=PREDICTIVE;  1=GENERATE  RCDF  VALUES] 
GENERATE  FULL  PRINTOUT  INDICATOR  [0=NO;  1=YES] 

START-UP  DYNAMICS  INDICATOR  [0=NO;  1=YES] 

TURBINE#  [0=INPUT  LIST;  3=TURBINE  3;  7=MIT  ACE;  8=FULL  ACE] 

If  ITURB=0,  insert  the  following  turbine  parameters  (one  per  line) . 


PIDP 

DESIRED  PRESSURE  RATIO 

ETADP 

ADIABATIC  EFFICIENCY 

GDP 

SPECIFIC  HEAT  RATIO 

AT 

EFFECIVE  CHOKED  AREA  (SQUARE  INCHES) 

IR 

ROTOR  MOMENT  OF  INERTIA 

(KG  M**2) 

DTIP 

TIP  DIAMETER  (INCHES) 

TGTWDP 

GAS  TO  WALL  TEMPERATURE 

RATIO 

REDP 

DESIGN  REYNOLDS  NUMBER 

DESIGN  OPERATING  SPEED  [RPM] 

EDDY-BRAKE  OPTION  [1-COMPUTED  BY  MODEL,  2-PREPROGRAMMED] 


If  IBEX=2,  insert  the  following  parameters  (one  per  line) . 

BO  INITIAL  EXCITATION  LEVEL 

DB_DT  SECOND  ORDER  EXCITATION  CONSTANT 

D2B_DT2  THIRD  ORDER  EXCITATION  CONSTANT 

SUPPLY  TANK  PRESSURE  (PSIA) 

SUPPLY  TANK  TEMPERATURE  (+K 1  -F) 

DUMP  TANK  PRESSURE  (PSIA) 

DUMP  TANK  TEMPERATURE  (+K I  -F) 

DESIRED  TEST  GAS  GAMMA 

If  IFRCDF=1/  input  a  single  line  with  the  following  four  values: 
PFILL(l)  PARTIAL  PRESSURE  OF  GAS  #1  (PSIA) 

TFILLK(l)  TEMPERATURE  OF  GAS  #1  (K) 

PFILL(2)  FINAL  FILL  PRESSURE  (PSIA) 

TFILLK(2)  FINAL  FILL  TEMPERATURE  (K) 
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Outline  of  Response  File  Entries  for  the  Facility  Model  Program  (continued) 

Parameter  Usage  in  Facility  Model  Program 

NTGCMP  NUMBER  OF  TEST  GAS  COMPONENTS  [NTGCMP  =  1  OR  2] 

ITGCMP  TEST  GAS  INDICATOR(S)  [1=AIR;  2=ARGON;  3=FREON12;  4=N2;  5=C02] 

P5PSI  OPENING  RESERVOIR  PRESSURE  (PSIA) 

P8PSI  CLOSING  RESERVOIR  PRESSURE  (PSIA) 

PVRPSI  VENT  RESERVOIR  PRESSURE  (PSIA) 

PVAOPSI  INITLAL  ACTUATOR  PRESSURE  (PSIA) 

NAGCMP  NUMBER  OF  ACTUATOR  GAS  COMPONENTS  [USUALLY  1] 

lAGCMP  ACTUATOR  GAS  INDICATOR  [  1=AIR;  2= ARGON;  3=FREON  12;  4=N2;  5=C02] 

TSIM  SIMULATION  RUN  TIME  (SECS) 

NPMS  NUMBER  OF  STEPS  PER  MILLISECOND 

TOPEN  VALVE  OPENING  TIME  (SECS) 

TCLOSE  VALVE  CLOSING  TIME  (SECS) 

TON  BRAKE  TURN-ON  TIME  (SECS) 

TOFF  BRAKE  TURN-OFF  TIME  (SECS) 

Hopefully,  this  document  provides  an  accurate  summary  of  the  input  data  required  to  success¬ 
fully  use  program /fld/z  O'-  It  is  important  to  note  that  when  the  program  is  run  in  the  mode  to 
generate  RCDF  values,  it  is  essential  that  a  file  of  the  form  "run*.rcdf"  (*  =  run  number)  be  pre¬ 
sent  in  the  current  working  directory.  No  such  restriction  is  imposed  when  the  program  is  run  in 
predictive  mode. 
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Technical  Memo 


To:  Dr.  C.  MacArthur 
From:  C.  Haldeman,  M.  Dunn 
Subject:  Main  Valve  Model  Update 
Date:  7-27-92 


Dr.  MacArthur- 

The  purpose  of  this  brief  note  is  to  bring  you  up-to-date  on  our  progress  with  the  main 
valve  model,  show  you  some  results,  outline  the  testing  intended  for  this  week,  and  discuss 
quickly  some  concerns  so  that  you  can  be  thinking  about  them. 

First  of  all  we  have  a  model  which  predicts  the  X  vs  T  of  the  main  valve  well,  and  it 
predicts  the  pressures  in  the  chambers  (although  not  quite  as  well).  This  model  was  developed  by 
examining  one  ran  of  a  test  matrix  in  detail  and  then  using  the  best  model  from  that  one  ran  to 
predict  other  runs,  varying  only  those  parameters  which  varied  from  ran  to  run.  Thus  things  like 
the  friction  and  grove  placement  would  remain  the  same  for  both  runs  and  only  the  valve  areas 


would  change.  The  p^cular  matrix  we  were  examining  was  test  matrix  3  consisting  of: 


Run 

PI 

Open 

Res. 

(psia) 

P2 

Act. 

Res. 

(psia) 

P3 

Close 

Res. 

(psia) 

Notes 

1 

331 

174 

426 

Both  Vent  valves  removed.  All  other  valves  set  to 
10.5  turns  open 

2 

335 

181 

443  n 

Same  as  Run  one  but  all  valves  wide  open 

3 

335 

180 

432 

Same  as  Run  2,  but  Vent  lines  removed 

The  real  purpose  of  this  test  was  to:  1)  examine  the  effect  of  the  smaller  reservoir  installed  on  the 
opening  side  and  2)  see  if  the  vent  lines  were  restricting  the  valve  motion. 

Run  1  was  used  as  the  design  ran,  with  all  the  parameters  being  set  to  achieve  "optimal" 
agreement,  and  then  this  model  was  used  to  predict  runs  2  and  3.  By  way  of  summary,  the  model 
has  been  changed  in  three  key  areas: 

Actuator  volumes: 

Original: 

V(6)=AO*S+VF6 

VC7)=AC*(SMAX-S)+VF7,  SMAX  =12 

New: 

V(6)=AO*S+VF6+VDEAD 
V(7)=AC*(SMAX-S-SP)+VF7+VDEAD 
Where  SMAX  =  internal  length  of  main  valve  piston  chamber  =  15.5 
AO=AC=  Cross  sectional  area  of  piston 
S  =  position  of  Actuator  piston 
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VF7,  VF6  =  the  volumes  of  the  respective  fill  lines 
VDEAD  =  Volume  of  the  head  plate 

Cross  Sectional 

Area  AO=AC  Opening  side 
Closing  side _  , _ (VDEAD) 


1  ^ 

1 

1 

f 

- —I - 

1 

1 

_ 1 - 

1 

1 

1 

1 

1 — 

1 

1 

1 

I 

! 

I 

1 

^  A  V 

Close  side  inlet  (VF7)  Open  side  inlet  (Vr6) 

The  key  difference  here  is  the  addition  of  the  volume  (9.6  in^)  under  the  head  which  was  not 
oiisinally  accounted  for .  Previously,  I  thought  that  Gerry  did  not  account  for  the  piston  thickness 
(3.5").  This  was  because  I  thought  he  was  using  SMAX  =  15.5;  however,  he  ran  his  models 
with  a  different  SMAX  ( =  12)  which  intrinsically  accounts  for  the  piston  thickness. 

Frictional  Coefficients: 

Originally  the  model  had  only  two  friction  values;  a  static  and  dynamic  friction.  I  have 
found  that  a  more  complicated  scheme  provides  better  agreement  This  is  a  piecewise  continuous 
system 

Friction 


CF| 

1 

CD2 

CD3 

CD4 

CDl 

1 

0  XI  X2  X3  10.125 

which  has  a  basis  in  the  fact  that  the  the  friction  is  reduced  dramatically  when  the  piston  rings  are 
over  the  grooved  area  of  the  cylinder. 

Actuator  Bypass  Area 
Original: 

If  (S<  SI)  Then  A67=AF67 
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If  (S>S1  and  S<S2)  Then  A67  =  AF67  +AV67 

If  (S>S2)  Then  A67=  AF67 

Sl=5,  S2=10.5,  AF67  =  Cross  over  area  (outside),  AV67  =  Area  of  grooves 
New: 

If  (S  <  SO  or  S  >SSW)  Then  A67=  AF67 

If  (S  >  SO  and  S<  S 1)  Then  A67=  AF67+  ALEAK 

If  (S>S1  and  S<S2)  Then  A67  =  AF67  +AV67 

If  (S>S2  and  S  <SSW )  Then  A67=  AF67  +  ALEAK 
Where  ALEAK  is  an  extra  leakage  which  stops  when  the  valve  stops  moving  at  either  end  and 
SSW  is  a  variable  (usually  fixed)  at  the  maximum  displacement  of  the  valve  (1 1.125  "). 

Results: 

The  main  criteria  for  a  good  model  were,  in  the  order  of  preference: 

1)  Match  of  X  vs  T 

2)  Match  of  chamber  pressures  after  valve  stops 

3)  Match  of  reservoir  pressures 

4)  Similarity  between  closing  and  opening  characteristics 

5)  Match  of  chamber  pressures  during  valve  motion 

Hgure  1  shows  the  X  vs  T  for  all  three  runs,  the  differences  in  closing  time  is  due  to  the  manual 
firing  of  the  main  valve.  Figure  2  shows  the  model  calculated  from  run  1  (called  model  96)  used  to 
predict  the  motion  of  run  2.  A  couple  of  important  points.  First  the  model  shows  the  compression 
of  the  springs  (which  we  do  not  measure),  thus  the  model  shows  a  greater  final  X  position. 
Secondly,  the  two  places  of  greatest  discrepancy  occur  at  the  same  relative  position  in  the  stroke, 
after  about  9  inches  of  movement  Finally,  it  is  important  to  realiK  that  after  the  valve  is  opened 
8.5"  it  is  fuUy  open,  any  more  distance  does  not  change  the  flow  area.  The  sources  of  these  . 
inconsistencies  in  the  model  are  shown  in  figure  3.  Both  on  the  opening  and  closing  side,  the 
model  over-predicts  the  dampening  pressure.  I  am  not  sure  why  this  occurs.  Figure  4  shows  the-, 
reservoir  pressures  which  agree  well. 

The  last  major  inconsistency  in  the  model  is  this  over-prediction  on  the  dampening  side  of 
the  piston  and  the  offset  in  the  dampening  side  (i.e.  closing  side  for  the  opening  stroke  and 
opening  side  for  the  closing  stroke)  after  the  piston  stops  moving.  Any  corrections  to  account  for 
this  make  the  X  vs  T  change. 

For  con^arison  purposes,  figures  5  and  6  show  the  output  from  the  original  math  model 
with  the  same  areas  and  volumes  that  I  use  in  my  models  but  with  the  contributions  of  ALEAK 
grouped  into  AF67  and  with  no  VDEAD.  One  can  see  that  the  original  model  under-predicts  both 
position  and  pressure  differences  greatly  during  the  valve  motion. 

V  Future  Work 
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Using  the  model  as  it  now  exists,  two  predictions  have  been  made.  Figure  7  shows  the 
first  run  which  will  be  nied  with  pressure  in  the  supply  tank.  This  test  corresponds  to  15  psig 
across  the  valve.  One  can  see  that  the  velocity  is  low  (.25  m/s)  as  the  valve  closes  (nominal  run 
conditions  according  to  Gerry's  model  were  about  .75  m/s).  Figure  8  shows  a  nominal  high  AP 
run.  One  important  point  to  note  is  the  vibrations  which  occur.  This  is  due  to  the  fact  that  it 
requires  so  much  initial  pressure  to  overcome  the  supply  tank  pressure  that  once  this  occurs,  the 
piston  goes  extremely  fast,  compressing  the  gas  in  the  dampening  side.  However,  we  do  know 
that  the  model  as  it  presently  exists,  over-predicts  these  pulses,  and  thus  these  vibrations  may  in 
fact  be  much  less.  Even  if  they  did  occur,  they  are  occurring  at  a  location  where  the  main  valve  is 
not  influencing  the  flow  area  (after  about  8.5").  And  finally  the  velocities  seem  to  be  reasonably 
low  during  the  final  stages  of  motion. 

With  this  in  mind,  the  goal  of  these  week  is  to  steadily  increase  supply  tank  pressure  while 
measuring  valve  performance.  At  this  point  we  have  reason  to  believe  that: 

1)  The  model  does  predict  pressures  and  position  fairly  well 

2)  The  possibility  for  damage  to  the  valve  during  this  ramp  up  process  is  minimal 

There  is  stiiU  work  which  needs  to  be  be  done  on  the  model.  The  two  areas  which  need  to 

be  addressed  are  the  pressure  inconstancies  mentioned  earlier  and  the  questions  involving  the 
proper  prediction  of  velocity.  The  model  uses  a  AX/AT  system  to  predict  velocity.  This  is  clearly 
an  approximation  to  the  real  velocity  which  is  8X/9T.  Presently  an  inconsistency  exists  between 
the  AX/AT  data  and  the  data  which  is  obtained  by  fitting  a  polynomial  through  X  vs.T  and  then 
analytically  differentiating  it  To  some  degree  this  problem,  while  annoying,  is  not  major  since  the 
impact  criteria  can  be  set  for  velocities  in  either  system.  Presently  our  design  criteria  is  not  to  hit 
the  0-ring  sealing  plate  any  harder  then  we  have  been,  and  in  no  case  to  hit  it  such  that  the  springs 
compress  entirely.  At  this  point  any  future  iir^rovements  in  the  model  can  be  relegated  to  a  future 
task,  when  other  more  critical  issues  have  been  resolved. 
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To:  C.  MacArthur  and  R.  Bergerly 
From:  C.  Haldeman 
Subject:  Math  Model  Revisions 
Date:  10-27-92 


I  have  become  convinced  that  during  valve  motion  there  is  something  physically  happening 
which  we  are  not  modeling  well.  This  comes  from  two  important  facts: 

1)  The  long  term  pressure  traces  (i.e  over  2-5  seconds)  agree  well  for  the  reservoir 
pressures,  and  the  patterns  are  similar  for  chambers  6  and  7  (although  the  offsets  are  off) 

2)  Any  attempt  to  correct  the  short-term  pressure  distribution  in  chambers  6  and  7  (i.e 
while  the  piston  is  moving)  degrades  our  relatively  good  agreement  with  the  reservoir  pressures. 


One  explanation  of  this  is  that  during  valve  motion  another  choking  condition  appears. 
This  would  explain  the  high  back  pressure  in  the  compression  chamber.  The  figure  below  shows 
where  I  believe  these  two  choke  conditions  occur. 


AV67 


Existing  chokes 


Proposed 

chokes 


There  are  two  ways  to  proceed.  The  most  complete  way  would  be  to  remodel  the  valve  system  as 
comprising  of  6  chambers  instead  of  the  existing  4.  However  that  could  become  quite 
complicated.  A  simpler,  quicker  solution  may  be  to  change  the  venting  areas  and 
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compression  volumes  when  the  piston  is  moving.  Thus,  on  the  opening  stoke  when  the  piston 
is  moving,  the  total  volume  of  chamber  7  would  only  be  chamber  7  and  not  the  line  volume  and 
VDEAD,  and  the  choke  area  would  be  different  from  A7V.  After  the  piston  stops  moving  the  total 
volume  increases  by  the  line  volume  and  VDEAD  and  the  choke  are  changes  to  A7V.  On  the 

closing  stroke  the  same  thing  would  happen  to  chamber  6. 

We  measure  the  chamber  pressures  (6  and  7)  outside  of  the  main  valve  itself,  and  thus, 
outside  of  the  new  proposed  choke  areas.  This  would  imply  that  our  measured  values  on  the 
compression  side  are  acmaHy  lower  than  what  would  be  in  the  chamber  during  this  choking 
condition.  However  this  could  explain  quite  a  bit  about  the  compression  spike.  Looking  at  figure 
9  on  the  last  main  valve  memo  I  sent  (10-6-92),  one  sees  that  the  model  predicts  a  large 
compression  spike  which  we  do  not  see.  K  there  was  a  secondary  choke  occurring  this  could  in 
effect  isolate  our  compression  pressure  transducer  from  measuring  this  spike.  This  scenario  would 
also  explain  the  strange  friction  factors  we  obtain  from  our  measurements  since  the  compression 

side  pressure  would  not  be  accurate.  Here  are  my  suggestions; 

1)  Use  timing  information,  rather  than  valve  position  to  set  when  the  choke  areas  are 

activated. 

2)  Assume  that  the  choke  starts  in  front  of  the  head  (so  VDEAD  is  added  to  the  line 
volume). 

3)  Make  all  these  parameters  controllable  from  cwh_config.dat 
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To:  Dr.  C  MacArthur,  C.Murawski,  and  J.  Finnegan 
From:  C.  Haldeman,  M.  Dunn 
Subject:  Main  Valve  Model  Update  II 
Date:  10-6-92 

The  purpose  of  this  brief  note  is  to  outline  the  results  of  the  main  valve  model  for  the  high 
pressure  runs  (Matrix  6),  describe  the  stams  of  the  model,  and  to  transfer  to  WPAFB  the  capability 
of  running  the  model.  The  present  changes  in  the  model  are  outlined  in  the  memo  dated  7-27-92  to 

Dr.  Mac.ALTthur. 

Nomenclature: 

The  model  results  conform  to  Dr.  Guenette's  ori^al  notation.  Figure  1  is  a  schematic  of 
what  the  different  areas  and  chambers  represenn  V5  is  the  open  reservoir  and  Vg  is  the  closing 
reservoir.  The  five  choke  areas  A7V,  A87,  .AF67,  A6V,  and  A56  are  physically  created  by  the  five 
needle  valves.  .AV67  represents  the  grooves  cut  into  the  piston  cylinder  housing.  The  following 
terms  are  used  in  the  model  and  are  sometimes  used  on  the  graphs: 

POS-  Valve  position,  sometimes  referred  to  as  X,  units  are  inches 
VEL  -  Valve  velocity,  unites  are  m/s 

Px-  Where  x  is  the  chamber  number,  refers  to  the  pressure  in  the  chamber,  units  are  psi 
Throuahout  the  write-up  there  are  several  terms  used  to  describe  the  main  valve  behavior  and 
pressure  fluctuations  in  chambers  6  and  7.  These  are  shown  on  figure  4  and  are  summarized 

below. 

Opening  Stroke-  The  operation  of  the  valve  going  from  closed  to  open  position 
Closing  Stroke-  The  operation  of  the  valve  going  fi'om  open  to  closed  position 
Initial  Pressure  Spike  -  The  first  pressure  pulse  in  the  acraating  side  of  the  piston  (P6  during 
the  opening  stroke  and  P7  during  the  closing  stroke) 

Compression  Pressure  Spike  -  The  first  pressure  pulse  in  the  compression  side  of  the  piston 
(P7  during  the  opening  stroke  and  P6  during  the  closing  stroke) 

Predicted  Performance  of  Main  Valve: 

Figures  2-6  show  the  predicted  performance  of  the  main  valve  (for  the  present  model 
configuration)  for  differential  pressures  of  50  psi,  80  psi  and  95  psi.  Figure  2  shows  the  position 
and  velocity  prediction  of  the  valve.  Hgure  3  shows  the  pressure  histories  of  the  opening  and 
closing  reservoirs.  Figure  4  shows  the  pressure  histories  of  volumes  6  and  7  (this  figure  also  has 
examples  of  the  pressure  spikes).  Figure  5  and  6  are  the  position  and  velocity  curves  for  the  80 
and  95  psi  cases  respectively. 
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Model  Summary: 

Figures  7-9  show  the  analysis  for  ±e  last  set  of  data  points  obtained  (the  55  psi  run,  Mamx 
6,  Run  6)  Fisure  7  shows  the  pressure  histories  for  the  open  and  close  reservoir,  while  figure  8 
shows  the  valve  position  and  figure  9  the  pressure  histories  inside  the  piston.  While  figure  7 
shows  good  agreement  between  the  recorded  reservoir  pressures  and  the  model  figure  8  shows 
that  the  piston  position  does  not  agree  as  weU.  OrigmaHy  the  model  initial  pressure  spike  (see 
fioure  9)  was  in  a  sisnificantly  different  time  position  relative  to  the  experimental  oaia.  This  was 
fo^d  to  be  a  timing  problem  which  was  easily  corrected.  While  the  model  stffl  shows  all  the  same 
characteristics  (i.e.  there  are  spikes  at  the  same  time  on  both  the  model  and  experimental  data,  the 
absolute  levels  differ  signincantly. 

The  initial  pressure  spike  is  a  very  strong  function  of  the  reservoir  choke  areas:  A87  or  A5 
(meaning  that  very  small  changes  in  the  area  make  the  model  either  over-predict  or  under-predict 
this  phenomenon).  .And  from  the  valve  position  histories,  it  is  clear  that  the  absolute  values  have 
little  influence  on  the  valve  motion.  It  is  the  compression  spike  that  is  causing  prooiems  for  the  ^ 
valve  motion.  Several  attempts  have  been  made  to  adjust  the  model  to  make  the  compression  spike 
more  compatible  with  experimental  results  (discussed  later),  but  none  were  sucw,.ssfai.  At  this 
point,  the  model  can  be  used  predictive  tool  to  estimate  ±e  valve  motion.  Since  we  have  shown 
that  the  system  is  over-damped,  the  chances  of  damaging  the  valve  have  been  reauceo  considerably 
and  it  is  our  recommendation  to  continue  with  the  high  pressure  testing. 

Model  .4nalysis: 

Several  attempts  were  made  to  make  the  compression  spike  more  closely  approximate  ±e 
experimental  data.  The  spike  could  be  reduced  by  increasing  the  leak^e  area  {ALEAK)  to  about  .1 
in2  (which  is  about  2.5  times  the  flow  area  of  the  gas  leaving  the  reservoirs).  However,  the 
‘  predicted  decaying  pressure  history  of  the  closing  acmator  (from  .8  to  2.1  seconds,  se«  figure  9) 
was  stiU  much  lower  than  the  experimental  one.  In  addition,  the  volume  of  the  closing  side  was 

decreased  dramatically  (with  no  influence).  ^  ■ - 

It  has  been  generally  concluded  that  the  size  of  the  reservoirs  and  the  opening  and  ventmg 

choke  areas  (A87,  A56,  A7V,  A6V)  set  the  macro  (0-5  secs)  pressure  decay  of  the  reservoirs  and 
the  timing  of  the  initial  pressure  spike.  The  slope  of  the  decay  (see  figure  7)  is  dictated  by  the 
internal  passages  in  the  piston  cylinder.  TTie  slope  of  the  decay  of  the  inner  chambers  (see  figure 

9)  is  dictated  by  the  bypass  ratio,  AF67. 

Since  the  valve  motion  is  dictated  by  the  difference  between  ?6  and  P?  in  the  first  second  of 
operation,  the  pressures  after  that  may  not  need  to  line  up.  However,  pressures  that  are  far  off  ^ 
when  the  valve  goes  to  close,  cause  much  larger  differences  between  prediction  and  experiment  lor 

the  closing  stroke. 
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Running  the  Valve  Model: 

It  is  clear,  that  the  best  way  to  improve  valve  performance  is  to  integrate  ±e  model  testing 
with  valve  testing.  The  most  inexpensive  and  expedient  way  to  do  this  is  to  have  the  people  at 
WP.^FB  do  the  integration.  The  model  is  stored  in  Haldeman’s  account  and  can  be  accessed 
simply  by  typing  "facility"  from  that  account  It  requires  two  input  files.  One  is  the  configuration 
file*"cwh_config.dat"  (a  copy  is  enclosed).  This  file  has  aU  the  information  about  the  main  valve 
such  as  the  choke  areas,  the  reservoir  volumes,  etc.  The  second  file  "facility.startup  contains  aH 
the  initial  conditions  of  the  facility  such  as  the  supply  tank  pressure,  opening  time  of  valve,  etc. 
Thus  if  the  valve  were  working  perfectly,  one  would  only  need  to  modify  this  last  file  to  run  any 

test  condition. 

The  model  has  three  types  of  output  There  is  the  traditional  output  (labeled  BD.OUT  and 
BD.PLT)  which  contain  aU  ±e  model  information.  BD.OUT  also  contains  the  initial  information 

stored  in  "cwh_config.dat"  and  "facility.semp".  BD.PLT  has  its  data  comma  delimiKd  Anew 

file  is  called  "valve.dat"  which  contains  information  specific  to  the  valve  (comma  delimited). 

These  three  files  are  ASCII  and  can  be  read  by  other  programs  or  computers  (like  a  Macmtosh). 

The  third  set  of  output  are  the  velocity,  pressure  and  positions  associated  with  the  main  valve  m 
DAS  format  These  are  called  P5_FM.3,  P6_FM.3,  etc.  These  files  can  be  used  as  input  to  the 
DRP,  so  it  is  easy  to  make  comparisons  between  the  model  and  the  experimental  data. 

There  are  three  Fortran  routines  which  help  in  data  processing.  One  is  called  dIvfmtF 
which  converts  from  DAS  to  ASCH  format  Tnis  program  can  be  used  to  take  several  DAS  files 
and  create  a  multiple  column  ASCH  file  (which  is  quite  large).  The  correct  usage  can  be 
determined  by  just  typing  the  program  name,  and  it  responds  with  the  appropriate  smtax.  Once 
this  ASCn  file  has  been  aeated,  another  program  called”timintF'  can  be  used  to  select  a  piece  of 
this  bigger  file  based  on  time  (i.e.  aH  data  between  1  and  2  seconds)  or  select  every  n^^  data  point 
The  output  can  then  usually  be  handled  by  commerciaUy  available  analysis  programs.  The  final 
proaram  is  called  "foraftF’  which  takes  as  its  input  the  two  LED  files  used  to  measure  the 
position  of  the  main  valve.  This  program  finds  the  comers  of  the  pulses  and  is  used  to  determine 
when  the  main  valve  reverses  direction.  The  output  is  also  ASCn  and  can  be  ported  to  other  data 
analysis  programs.  This  has  been  used  as  input  to  the  Macintosh  which  actuaHy  develops  the 

position  of  the  valve  over  time. 

Future  Work: 

Preliminary  predictions  suggest  that  if  ah  needle  valves  were  eliminated,  then  the  main 
valve  may  move  smoother  at  higher  pressures  (i.e.  with  less  vibration).  Choke  areas  are 
approximately  .06  in^  without  the  needle  valves,  .04  in2  with  the  needle  valves  wide  open,  and 
.025  in2  with  the  needle  valves  ten  turns  open.  All  values  shown  in  "cwh_config.dat"  are  the  best 

ones  known  to  date. 
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We  believe  that  the  present  configuration  will  work  at  ail  pressure  values,  although  we  also 
beHeve  that  the  system  could  be  optimized  by  tweaking  some  of  the  valves  on  the  closmg  ade. 
Rsure  8  shows  that  the  main  valve  covers  about  90%  (9  inches)  of  its  closing  stroke,  reladvely 
quickly,  while  the  last  10%  (1")  takes  150  ms.  This  should  be  capable  of  being  improved  upon. 

It  is  clear  that  something  is  happening  in  the  con^ssion  chamber  which  causes  the  ^ 
pressure  to  be  quite  high,  to  that  predicted  by  the  model,  and  that  this  is  not  a  function  of 

the  vent  or  inlet  choke  areas  (otherwise  the  reservoir  pressure  histories  would  not  a^).  ^ 
Understanding  this  may  take  some  extra  experimentation,  and  possibly  some  derivation  of  internal 
properties  based  on  measured  values.  Presently,  this  can  best  be  done  by  the  personnel  at  WPAFB 
who  have  access  to  the  main  valve.  We  are,  of  course,  ready  and  willing  to  help  n  any  difficulties 

should  arise. 
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cwh__config.dat 
*/  Facility  Configuraiion  Data  File 
*/  FDl  has  Facility  Volume  Data,  ft**3 
W- SUPPLY  TANK 
-  TURBINE  INLET 
*/V2  -  TURBINE  EXHAUST 
*A^3-DSV  INLET 
*A^4- DUMP  TANK 
SFDl 

V0=  3200., 

Vl=  21., 

V2=  2., 

V3=  250., 

V4=  6400. 

SEND 

*/FD2  has  Facility  Flow  Path  Orifice,  in**2 

*/ADSO  -  EXIT  PIPES  TO  DUMP,  A34: 2-20"DIA,CD=.62 

SFD2 

ADSO=  390. 

SESD 

* 

*/VDl  has  Valve  Geometry  variables 

*yiSS  =1  FOR  SLIDING  SLEEVE  (ATAR),  =0  FOR  COWCATCHER  (MTT) 
*/RO  -  VALVE  ANNULUS,  OUTER  RADIUS  (INCHES) 

*/RI  -  VALVE  ANNULUS,  INNER  RADIUS  (INCHES) 

*/RM  -  MAIN  SEAL  RADIUS  (INCHES) 

*/XM  -  MAIN  SEAL  BREAKAWAY  (INCHES) 

*/KS  -  BACK  SEAL  RADIUS  (INCHES) 

*/XS  -  BACK  SEAL  BREAKAWAY  (INCHES) 

*yTHETA  -  ENTRANCE  ANGLE  FROM  VERTICAL  QDEGREES) 

*/MASS  -  TOTAL  MOVING  MASS  (KG) 

SVDl 
ISS=  1, 

RO=  17.25, 

RI=  7.75, 
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RM=  18.00, 

XM=  0.058, 

RS=  17.30, 

XS=  0.041, 

THETA  =  45.0, 

MASS  =  277.0 
SEND 

* 

*  VD2  has  Impact  Spring  Data 

*/XKSC  -  ZERO  FORCE  POSITION,  MAIN  VALVE  CLOSING  SPRING  (IN) 

*/KSC  -  SPRING  CONSTANT  (LBF/IN),  CLOSING  SPRING 
*/XKSO  -  ZERO  FORCE  POSITION,  VALVE  OPENING  SPRING  (IN) 

*/KSO  -  SPRING  CONSTAITT  (LBF/IN),  CLOSING  SPRING 
*/CSF  -  COEFFICIENT  OF  STATIC  FRICTION 

* /  COEFFICIENTS  OF  DYNAMIC  FRICTION  FOR  SLIDER  RANGES  SHOWN: 

*/CDl-0TOXl 

*/CD2-Xl  TOX2 

*/CD3-X2TOX3 

*/CD4-X3TOENT) 

*/Xl  -  FIRST  SLIDER  POSITION  FOR  DYNAMIC  FRICTION  COEFFICIENTS.  (IN) 
*yX2  -  SECOND  SLIDER  POSITION  FOR  DYNAMIC  FRICTION  COEFFICIENTS  (IN) 
*/X3  -  THIRD  SLIDER  POSITION  FOR  DYNAMIC  FRICTION  COEFFICIENTS  (IN) 

$VD2 


XKSC= 

=  0., 

KSC= 

8000., 

XKSO: 

=  10.125, 

KSO= 

5000., 

CSF= 

0.40, 

CD1= 

0.27, 

CD2= 

0.4, 

CD3= 

0.4, 

CD4= 

0.27, 

11 

2.0, 

X2= 

5.0, 

X3= 

8.0, 

SEND 
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*  VD3  has  Actuator  Cylinder  Geometry 
*/DB  -  BORE  DIAMETER  (EN.) 

*/DR  -  ROD  DIAMETER  (IN.) 

*/NRODS  -  =2  FOR  ATAR,  =1  FOR  MU 

*/SMAX  -  MAXIMUM  ACTUATOR  STROKE  (INCHES) 

*/SP  -  PISTON  THICKNESS  (INCHES) 

*/S0  -  STROKE  OFFSET,  S  @  X=0  (INCHES) 

*/S  1  -  OFFSET  TO  START  OF  TRANSFER  PORT  (IN.) 
*/S2  -  END  OF  TRANSFER  PORT  (INCHES) 

*/AV67  -  TRANSFER  PORT  AREA  (IN**2) 

*/AF67  -  ACTUATOR  FIXED  BLEED  AREA  (IN**2) 
*/AT  .F.AK .  LEAKAGE  AROUND  PISTON  (IN**2) 
*/SSW  -SWUCHPOSmON(IN) 

SVD3 

DB=  8.0, 

DR=  3.5, 

NRODS=  2, 

SMAX=  15.5, 

SP=  3.5, 

S0=  1.0, 

Sl=  4.75, 

S2=  7.75, 

AV67=  2.0, 

AF67=  0.030, 

ALEAK=  0.015 
SSW=  11.125, 

SEND 


*  VD4  has  Actuator  Gas  Supply  Geometry 

*/W5  -  OPENING  GAS  SUPPLY  VOLUME  (M**3) 

*/VF6  -  OPENING  LINE  VOLUME  (IN**3) 

*/VF7  -  CLOSING  LINE  VOLUME  (IN**3) 

*m  -  CLOSING  GAS  SUPPLY  VOLUME  (M**3) 
*/IVNT  -  VENT  TO  DUMP=0,  VENT  TO  ROOM=l 
*/A056  -  OPENING  GAS  SUPPLY  ORIHCE  (IN**2) 
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*/A07V  -  CLOSING  GAS  SUPPLY  ORIHCE  (IN**2) 

*/AC87  -  CYLINDER  CLOSING  END  VENT  ORIFICE  (IN**2) 
*/AC6V  -  CYLINDER  OPENING  END  VENT  ORMCE  (IN**2) 
*/VDEAD  -  DEAD  VOLUME  IN  CYLINDER  HEAD  (IN**3) 
SVD4 

V5=  2472.0E-6, 

VF6=  332, 

VT7=  32.0, 

V8=  4206.0E-6, 

IVNT=  1, 

A056=  0.0360, 

A07V=  0.060, 

AC87=  0.042, 

AC6V=  0.064, 

\T5EAD=  9.6, 

SEND 

* 

*  3RK  has  BRAKE  CONSTANTS 

*/KB  -  BR.AKE  CONSTANT  (W-s**2  /  Tesia**2) 

*/NB  -  brake  speed  @  maximum  torque  (ipm) 

*/MDRUM  -  BRAKE  DRUM  MASS  (KG) 

*/CPD  -  BRAKE  DRUM  SPECIRC  HEAT  (J/KG-K) 

*/TAUB  -  BEAKE  TIME  CONSTANT  (SECONDS) 

SBRK 
KB=  80.829, 

N'B=  7461, 

MDRUM=  33.36, 

CPD=  427.0, 

TAUB=  0.01 
SEND 
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faciiity.setup 

0  /Generate  RCDF  values  ,l-yes,.0-no/ 

1  /Generate  full  printout,  1-yes,  2-no/ 

1  /Stan-uo  dynamics?  1-yes,  0-no/ 

3  /Turbine  number/ 

5992  /Operating  speed/ 

1  /Eddy-brake  option  1-coii^uted  by  model,  2-preprogrammed/ 
95.0  /supply  tank  pressure/ 

300.0  /supply  tank  temperature  K/ 

1.40  /Dump  tank  pressure/ 

300  /Dump  tank  temp./ 

1.4  /Gamma/ 

1  /#  of  test  gas  components/ 

1  /I  -  air,  4  -  nitrogen/ 

500.0  /  open  resevoir  Press/ 

600.0  /close  reservoir  pressure/ 

14.30  /  vent  pressure/ 

14.30  /  Actuation  initial  pressure/ 

1  /#  of  actuator  gas  components/ 

4  /gas  type  l-air,  2-argon,  3-R12, 4-N2, 5-C02/ 

5  /run  time/ 

10  /number  of  timesteps  per  ms/ 

.097  /open  time/ 

228  /close  time/ 

.089  /Brake  rum-on  time/ 

2.5  /Brake  tum-off  time/ 
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